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We present a concrete and explicit construction of a new scalar constraint operator for 
loop quantum gravity. The operator is defined on the recently introduced space of partially 
diffeomorphism invariant states, and this space is preserved by the action of the operator. 
To define the Euclidean part of the scalar constraint operator, we propose a specific regular¬ 
ization based on the idea of so-called ’’special” loops. The Lorentzian part of the quantum 
scalar constraint is merely the curvature operator that has been introduced in an earlier 
work. Due to the properties of the special loops assignment, the adjoint operator of the non- 
symmetric constraint operator is densely defined on the partially diffeomorphism invariant 
Hilbert space. This fact opens up the possibility of defining a symmetric scalar constraint 
operator as a suitable combination of the original operator and its adjoint. We also show 
that the algebra of the scalar constraint operators is anomaly free, and describe the structure 
of the kernel of these operators on a general level. 


I. INTRODUCTION 


The canonical quantization o^eneral relativity has come a long way since the formulation of 
the Ashtekar-Barbero variables [l|, 1^. As a generally covariant theorv general relativity has its 
dynamics encoded in constraints. Loop quantum gravity (LQG) that is the incarnation 

of the mentioned quantization program, succeeded in defining a Hilbert space of kinematical 
quantum states, and implementing and solving the Gauss constraints, which encode the SU{2) 
gauge invariance, and the spatial diffeomorphism constraints [^. The scalar constraints are 
technically more involved because of their complicated expression in terms of the Ashtekar- 
Barbero canonical variables. 

The first rigorous proposal of a scalar constraint operator was introduced by T. Thiemann in 
[^, based on some concepts discovered by C. Rovelli and L. Smolin in 11]. The construction in¬ 
volves the volume operator 0 and uses a mathematical artifact to suppress the non-polynomial 
character of the constraints in terms of LQG variables. As a result, the constraint operator is 
gauge invariant and anomaly-free. This operator acts on the Hilbert space of diffeomorphism 
invariant states, but does not preserve this space due to the presence of the lapse function in 
the operator. 

Recently, a new Hilbert space of partially diffeomorphism invariant states was introduced 
14l |. In that article, it was shown that upon some changes in the Thiemann’s regularization of the 
scalar constraints, the resulting quantum operator would preserve J^^x- Moreover, the operator 
would still be anomaly free and there would be possibilities to define a symmetric constraint 
operator, making discussions of self-adjointness extensions and spectral analysis more accessible. 

In the present article, we explicitly implement the scalar constraints for LQG verifyiM the 
criteria discussed in 0 - We base our construction on ideas and concepts introduced in [§, 0 , 
0 to deal with the Euclidean part of the constraint, and the use of the curvature operator 
introduced in [lfil | to dehne the Lorentzian part. The article is organized as follows. In section HIl 
we briefly review the classical Ashtekar formulation of general relativity. In section Hill we review 
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the Hilbert space of LQG, the implementation of SU{2) gauge invariance and the construction 
of the partially diffeomorphism invariant Hilbert space section IIVI we present the 

regularization of the classical scalar constraint allowing us to define a non-symmetric scalar 
constraint operator and its adjoint, both densely defined. We discuss the quantum algebra, the 
possibility of defining a symmetric constraint operator, then the solutions of the quantum scalar 
constraints; We close in section El with some comments and outlooks to future developments. 


II. CLASSICAL THEORY IN ASHTEKAR VARIABLES 


The 3 + 1 Hamiltonian formulation of general relativity, written in terms of the Ashtekar- 
Barbero variables 0,3 {A\,Ef) (the spatial index a and the su{2) index i take the values 
1, 2,3), manifests as a constrained SU{2) gauge theory. The spatial variable and its conjugate 
momentum E^, the densitized triad, verify the canonical relations 


{Alix),E^^iy)} = kl35i5]8{x,y) 
{A*(x),A^,(y)} =0 = {E-{x),E){y)} 


( 11 . 1 ) 

( 11 . 2 ) 


where k = SvrG and /3 is the Immirzi parameter. 

The constraints obtained in this formulation consist of the Gauss constraints Gi{x) (gauge 
constraints), spatial diffeomorphism (vector) constraints Ca{x) and scalar constraints C{x). 
They are first class constraints and can be expressed as follows: 


Gi{x) = (daE^{x) + Ai{x)El{x 
Ca{x) = ^Fl,{x)E\{x) - A* (x)^,(x). 


(H.3) 


C{x) = 


1 


2kl3^ 


eijkEf{x)E^j{x)E^,{x) 
V^l detS(x)| 


+ (l - v^l detS(x)| R{x) j, 


where s = 1 in the case of spacetime with Euclidean signature and s = — 1 in the case of 
Lorentzian signature, the curvature of the connection AJj, and R is the Ricci scalar of the 
metric tensor Qat on the 3-dimensional manifold T, (the relation between Qab and the variable Ef 
is given by = EfE\/\ detE|). This form of the scalar constraints was proposed by Domagala 


19|, and was used also in our recent paper 1^. It is an alternative to Thiemann’s form of the 
scalar constraint 0] used in [l3 |. 

Imposing the constraints pi.3p is equivalent to imposing their smeared versions 
e?(A) = J cfxA\x)giix) , CiN) = J d^xN’^{x)Ca{x) , C{N) = j d^xN{x)C{x), (H.4) 

S S E 

where A(x) = TiA^{x) is an arbitrary su{2) valued smearing function, while iV“(x) and N{x) 
are arbitrary real valued smearing functions called the shift and lapse respectively. 


The constraints algebra reads 

{^(A),^(A')} = a([A,A']), 
{a(A),C(iV)} = -^(/:^A), 
{a(A),C(iV)} = 0, 


{CiM),CiN)} = CiC^N), 
{C{M),C{N)} = C{£^N) , 


(11.5) 


{C(M), C(V)} = C{q'^\NM^b - NM^b]) + G{S{A, E)), 


where ^(A, E) is a certain function of the phase space variables, whose explicit expression is not 
relevant for this work but can be found in 0. 
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The quantization program of LQG is a canonical quantization following Dirac’s procedure. 
Namely, the phase space variables are quantized and a Hilbert space of functionals of the con¬ 
figuration variable A is constructed, then classical functions on the phase space are promoted to 
quantum operators and the constraints are imposed on the quantum level as operators equations 
in order to determine the physical Hilbert space. In the following section, we briefly present the 
construction of the Hilbert space in loop quantum gravity along with the implementation of the 
Gauss and spatial diffeomorphism constraints. 


III. LOOP QUANTUM GRAVITY: KINEMATICS 

Loop quantum gravity is an attempt to built a background independent quantum theory of 
gravity, therefore there is no reference to any background metric in defining the classical algebra 
to be quantized. Also, since the Poisson brackets are singular, we need to introduce 

smeared variables, holonomies and fluxes (defined below), obtained by integration of A and E 
respectively over appropriate submanifolds of S. 


A. Kinematical Hilbert space 


The kinematical space in LQG is defined as the space of cylindrical functions of the variable A, 
i.e., complex valued functions depending on the su(2)-valued differential 1-form A = A^Ti^dx^", 
where r* € su{2) is a basis of sm(2 ), through finitely many parallel transports (holonomies) 


he[A] 



(III.l) 


where e is an oriented finite curve (edge) in S. Then a kinematical quantum state T has the 
form 


T[A] = V^(he,[A],...,heJA]) (HL2) 

with a function ■0 : SU(2)"' —>■ C. The set 7 := {ei,..., e^} is called the graph of T. 

The space of all cylindrical functions with a graph 7 is denoted by Gyl..^ and the space of 
all cylindrical functions by Cyl. The kinematical Hilbert space of LQG, 'Hkinj is defined as the 
completion of Cyl with respect to the norm defined by a natural scalar product @] 

Hkin := ^ (111.3) 

While a connection operator “A” is not defined, every cylindrical function T also defines a 
multiplication operator 

(T^)T')[A] = T[A]T'[A]. (HL4) 

The derivative operator is the quantum flux operator, obtained by quantization of the flux 
corresponding to E, 

Ps,i ■= f ^dx’’ Adx'^eabcC{x)Ef{x), (HL5) 

'j s 

through an oriented 2-dimensional surface S' C S. Here ^ : S —)• su(2) is a (generalized) smearing 
function that may involve parallel transports depending on A. The flux operator corresponding 
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to the classical variable pil.5p is then 

(III. 6 ) 

xeS e 

where e runs through the germs^ beginning at x, and Ks{e) = — 1 , 0,1 depending on whether 
e goes down, along, or, respectively, up the surface S. The operator Jx,e,i is assigned to a pair 
(x, e). Its action on the function T € Cyl defined in (IIII.2h . with ei belonging to the germ e, is 
given by 


7 

de 


% he2 , •••, ^en,)) 


e=0 


(III.7) 


B. Gauss & spatial diffeomorphism constraints 

In order to complete the quantization program, it is necessary to implement the constraints 
pi.4p and solve them. The Gauss constraint operator can be easily defined in terms of fluxes, 
and its kernel is identified with the space of gauge invariant cylindrical functions 

/(^) = + 9 ~^dg), for every 5 € (7^(2, SU(2)). (III. 8 ) 

We denote their algebra (a subalgebra of Cyl) by Cyl*^, and the corresponding Hilbert space 
C A dense subspace of is spanned by the spin network functions. A spin network 

function is defined by a graph 7 with half integer (non zero) spins assigned to the edges, and 
SU{2) invariant tensors (intertwiners) assigned to the vertices. Then the space of all gauge 
invariant states can be written as the orthogonal sum 

= 0^? (ni.9) 

7 

where 7 ranges over all the classes of graphs^, and is the Hilbert space defined as the 
completion of the space Cyl^ spanned by the spin-network functions of graph 7 .^ 

Let us now turn to the vector constraint. Due to the absence of a well defined operator 
corresponding to the spatial diffeomorphism constraint functional, the construction of a space 
of diffeomorphism invariant states is achieved through a diffeomorphism averaging procedure 
j^. The elements of each of the sub-spaces are averaged with respect to all the smooth 
diffeomorphisms Diff°°(S) which map 7 into analytic graphs. Recall that given a diffeomorphism 
/ : S —)• S, its induced action on a cylindrical function pH.21 1 is 

{Uf^)[A] = V^(h^(,,)[A],...,h^(,„)[A]) . 

But since Diff°°(S) is a non-compact set and we do not know any probability measure on it, 
we have to define the averaging in CyT, the algebraic dual to Cyl. The resulting space is a 
Hilbert space of diffeomorphism invariant states, denoted with a scalar product naturally 

inherited from the scalar product on 


^ A germ beginning at a point x is the set of curves overlapping on a connected initial segment containing x. 

^ Two graphs 7 and 7' belong to the same class if 7' can be obtained from 7 by a sequence of the following 
moves: splitting of an edge, connecting two edges, changing a orientation of an edge. 

® An important subtlety is, that given a graph 7, we define spin-network functions by non-trivial representations 
of SU( 2 ) assigned to the edges of 7. In general, Cyl^ contains also spin-network functions defined by a graph 
7” obtained from 7 by removing one of the edges. 
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However, we know that a quantum operator corresponding to the scalar constraint C{N) in 
pi.4p would not preserve the Hilbert space because of the presence of the lapse function 

N. In other words, an operator C{N) is not invariant under spatial diffeomorphisms. This fact 
raises serious difficulties in the treatment of relevant questions such as self-adjointness, spectral 
resolution and anomaly-freeness of the constraints algebra. 

A solution to this issue was suggested recently in 1 ^ . It consists of introducing an interme¬ 
diate space, the vertex Hilbert space 'b’he idea is to construct from elements of the Hilbert 

space partial solutions to the vector constraints, by averaging the elements of each of the 
sub-spaces T-L^ with respect to all the smooth diffeomorphisms Diff°°(S)vert( 7 ) which act trivially 
in the set of vertices Vert( 7 ). Denote by TDiff°°(S).y the subset of Diff°°(S) which consists of all 
diffeomorphisms / such that /(y) = 7 and Uf acts trivially in T-L^, and by Diff“(S)vert( 7 ) tbe set 
of elements of Diff°°(S)vert( 7 ) which preserve the analyticity of 7 . The set of the transformations 
"Hkin induced by Diff^(S)vert( 7 ) can be identified with 


D., := Diff^(S)vert(7)/TDiff(S).,. 


(HI.IO) 


The averaging is defined in Cyl* through a rigging map 

r? : Cyl^ 

1^7) 

where N-y is the number of elements of D.y which preserve the graph 7 . 

The resulting //(T-j,) is a well defined linear functional on Cyl^. One then extends 7 ( 4 '-),) by 
linearity to the algebraic orthogonal sum (IIII.Qh . obtaining a map 

n-.ng^^CyW (HI.12) 

The vertex Hilbert space W^x then defined as the completion 

^vtx := = 05g, (HI.13) 

[7] 

under the norm induced by the natural scalar product 

= 7(4^) (7(4^0) • (ni.i 4 ) 

Each state in is invariant under the action of elements in Diff^(S)vert( 7 )- this sense, 
those states are partial solutions to the quantum vector constraint. They can become full 
solutions of the quantum vector constraint by a similar averaging with respect to the remaining 
diffeomorphisms Diff(S)/Diff(S)vert( 7 )) forming the space 


C Cyr (HI.ll) 

= jr Y. 


IV. LOOP QUANTUM GRAVITY: DYNAMICS 


The quantization of the scalar constraint we propose is carried out by treating separately 
the two terms of the constraint C{N), expressed in equation (lIV.ip below. The first term of 


C{N) (see (jll.Sp i. is quantized using the loop prescription introduced in [1^ to regularize the 


curvature of the Ashtekar connection, and Thiemann’s trick to remove the non-polynomial 
dependence on the canonical variables, caused by the presence of the factor 1/\J\ det £'(x)|. 
This is a special case of quantization of this term proposed in (l3 |. The new element is a 
specific, explicit proposal for the regulator. The second term of C{N) was already regularized 
and promoted to a quantum operator, the curvature operator (l 6 | . and we will go briefly through 
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the details of its construction below. Using that operator in our definition of the quantum scalar 
constraint is a true departure from the paper |l4l ]. 


A. Regularization of the scalar constraint 

The starting point is the expression 

1 /■ „ / eijkEf {x)E^Ax)FK{x) . ! - \ 

C{E) = ^l^d^x Nix) ( ^|detU(x) f^ ^ ^ V|detii;(x)| ii(x) j , (IV.l) 

In the case of s = 1 (the space-”time” signature + + ++), the choice /? = ±1 kills the second 
term (and corresponds to the original self-dual Ashtekar variables). For that reason we call the 
first term the ‘Euclidean’ part, and we call the second term the ‘Lorenzian’ part. 


1 . Euclidean part 

We first consider the Euclidean part. To express it in a non-singular form, we use Thiemann’s 
trick, which consists of using the identity 

(x)EJ(x) 2 u I, 

^ = -e“^"{A^(x), V}, (IV.2) 

3tE(x)| fc ^ 

where V is the volume of S, 


U := y d^x^\detE{x)\. (IV.3) 

s 

The Euclidean part C^{N) then takes the form 

'-=^1 d^^Nix)e-^-E^,ix){A'^,ix),V}, (IV.4) 

E 

The expression (IIV.4p is regularized via approximation of the integral by a Riemannian sum 
over a partition with e being a parameter characterizing the size of the cells A in by 

replacing A(x) with values of A at a point xa chosen in each cell A, and replacing the connection 

coefficients with parallel transports along open curves s/(A) and the curvature coefficients by 
the holonomies along loops a/j(A) 

Cl-W = E (IV,5) 

where h^') is the holonomy in a chosen SU{2) representation I and Wi = i-\/lil + 1)(2/ -|- 1) is a 
normalization factor^, the curves s/(A) and loops a/i^(A) are assigned to each cell A such that 
this functional converges to C^iN) in the limit e —>■ 0. Below we propose a specific assignment. 
But before introducing it in detail, we will remind another important element of the procedure. 



The representation I is left arbitrary in our construction. In representation I, we choose a basis = 1 , 2 , 3 ) 

of su{ 2 ), satisfying 

Tr(r«)=0, Tr(r«r«) = ^5i,. 
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The first, intermediate, step of the quantization is to dehne in a partition dependent 

quantum operator C^e{N). This operator will not have a limit when e —>■ 0. Still, by duality 
we want to obtain a well defined operator on that carries the diffeomorphism covariance 
property of the classical constraint. To accomplish that, we need to adapt our regulator to 
each graph 7 and the corresponding subspace Hy independently. We propose in this paper the 
following prescription: 

• is a triangulation, i.e. each cell A is a tetrahedron; 

• each tetrahedron A has at most one node of the graph 7 as one of its vertices; 

• each node v of the graph 7 coincides with a vertex of a tetrahedron A^ and xa„ = v; 

• if u is a node of 7 , then 

— u is a vertex of Uy tetrahedra A* saturating the neighborhood of v (i.e the tetrahedra 
meet at v and compose a closed neighborhood centered at u); 

— the edges of the tetrahedra A^ saturating the neighborhood of v do not overlap with 

the edges of 7 meeting at v, except for one tetrahedron, which we call The 

tetrahedron A^*^^ is adapted to one chosen ordered triple of edges (e/, ej, ck) 
meeting at v, i.e., the edges (s/, sj, sk) of A^*^^ meeting at v are segments of the 
edges (e/, ej, ex) of the graph 7 but do not coincide with them; 

— to the ordered triple of edges (s/, sj, sk) meeting at v there are assigned three 

loops oriented according to the order of the triple (s/, sj, sk)] 

— A loop verifies the following conditions: 

i. is an analytic curve; 

ii. lies in a surface defined through a canonical choice of coordinates adapted 
to the edges (s/, sj, sk) and does not intersect the graph® 7 at any point 
except at u; 

hi. is tangent to the two edges e/ and ej of the graph 7 at the vertex v up to 
orders kj + l and kj + 1 respectively, where ki{> 0 ) and kj{> 0 ) are respectively 
the orders of tangentiality of e/ and ej at the node®; 
iv. Denote by sij the edge of A^"^^ that links the edges (s/, sj) to form a triangle 
of the the tetrahedron The shape of the loop marries the shape of 

the triangle (s/, sj, s/j) as good as possible; 

This prescription for the adapted partition is twofold: The first part, which contains all the 
requirements except the conditions on the loops, coincides with some of the requirements on the 
partition in Thiemann’s approach to regularize the scalar constraint j^. In addition, in the 
number Uy is set to be equal to 8 for any node v of the graph thanks to a specihc procedure to 
construct the saturating structure around v. We could adopt the same procedure to fix Uy but 
it is a priori possible to keep it as a free parameter that is the same for all vertices, hence we 
drop the v label in rest of the article. 

The second part of the above prescription is about the conditions on the loop structure. We 
use a prescription, first introduced in [1^, different from the one in Thiemann’s construction in 
which the loop coincides with the triangle (s/, sj, sjj) of ■ The whole prescription 


® We do not show the construction of those coordinates nor the rooting procedure for the loop in this article, but 
we direct the reader to or for the details. 

® The order of tangentiality of an edge ei incident at a node v is the highest order of tangentiality of the edge ej 
with the remaining edges incident at v (see [T^l. 



is diffeomorphism invariant and it makes a loop assigned to a pair of edges unique up to diffeo- 
morphisms. As we will see later, the conditions on the loops also allow to introduce a densely 
dehned adjoint operator of the non-symmetric scalar constraint operator^, thereby providing a 
way to define a symmetric constraint operator (the key condition is that as in [3, 14] the loops 
do not overlap the given graph). In the rest of the article we refer to those loops as special loops. 

Having the adapted partition, we straightforwardly quantize the expression in (IIV.5P by re¬ 


placing the Poisson bracket of and V with 1/ih times the commutator of the corresponding 

operators, taking for V the internally regularized volume operator of 


V := l^ 

P 


Efe 


= 






8-3! 




x,ej,iJx,ej,jJx,eK,k 


I,J,K 


(IV. 6 ) 


where Ip is the Planck length, kq an overall averaging constant, ej runs through the set of germs 
starting at the point x, and e{ei,ej,eK) = sgn[det(e/, ej, ej^)]. Considering a gauge invariant 
state 'h.y with a graph 7 , the resulting operator acts as 


C§4N)^^ := Y, Ci{N) 

Ae'T'^ 

1 


^ ‘ Ae”?" vgAD'y 


.(0 


(IV.7) 


At this stage, the operator defined in pV.7p still depends on the triangulation The 
dependence on the triangulation is removed in three steps: 


a) Denote by R{v) the closed region formed by the n tetrahedra of saturating a 

vertex v. Here (...) contains the labels of the edges intersecting at v and defining a specific 
tetrahedron. Classically, as we take the limit e ^ 0 in the sense of refining the adapted 
triangulation to another adapted triangulation such that e' < e, we have 

(IV. 8 ) 

R{v) AiJK 



the label UK refers to one tetrahedron of R{v). In other words, the integral over R{v) 
converges to n times the integral over any tetrahedron of R{v) as we take the limit e —>■ 0. 
For the operator in pV.7p . this translates as 


ClPiV)^, := 


n 


E E 

\IJKci 






V 


(IV.9) 


b) A triangulation selects at each node u of a graph 7 a unique triple of edges (ej, ej, ex) 
meeting at v. In order to remove this selection from the operator, it is enough to average 
at each node v over the classes of triangulations that select different triples meeting at v. 
Therefore the operator would contain contributions from all possible triples meeting at 


^ In case of Thiemann’s construction, the adjoint operator of the non-symmetric scalar constraint operator is not 
densely defined. 
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the same node and we obtain 




n N{v) 

ihk‘W^ ^ E{v)' 


JJK 


Tr U 


,(0 




(0 


a/j(A) sk{A)L sk(A) 






dG7 


(IV. 10) 


where now the UK run through all triples of edges of the graph 7 meeting at the node v, 
and E{v) is the number of unordered triples of edges meeting at v (hence E{v) depends 
only on the graph 7 ). Notice that due to the presence of the volume operator in its 
expression, (7®^, annihilates two-valent nodes and nodes which have degenerate differential 
graph structure. Therefore the action of the operator Cf{N) on a gauge invariant state 
is always finite and it also preserves the gauge invariant space. 


c) The only dependence left on the triangulation is in e. We then need to take the limit 
e ^ 0. As we have mentioned above, in this limit Cf{N) does not converge to any well 
defined operator in the space The way around this problem is to first pass the 

operator Cf{N) to the space by duality, then take the limit 
is ensured and the final operator is then defined as 


14l |. The convergence 


C^(N) := lim 
£->■0 


Cf(iV) 


(IV.ll) 


acting in the space of gauge and partially diffeomorphism invariant states "H^x- 

The operator C^{N) is densely defined on the space 'H^x; it contains the span of partially 
diffeomorphism invariant spin networks space 7 (^), and graph changing as it removes special 
loops at the nodes®. It maps its domain C W^x to a subset of W^x therefore preserves 
the gauge and partial diffeomorphism invariance. 


2. Lorentzian part 

Now let us turn to the Lorentzian part of (IIV.1|) . namely 

l-s/32 


U{N) = 




./AiVWV|detE(x)|fl(x). 


(IV.12) 


The quantization of this classical functional was already carried out in 16(]. The regularization 
is external and based on the Regge approximation 17| of the 3d Einstein-Hilbert action. On 
a gauge invariant state the non-symmetric operator® corresponding to the Lorentzian part 
cicts clS 




l-s/32 

8k/3^ 




'ei,ej 


vGj 

= :^iV(u)a%, 


lU 


(IV.13) 


£>€7 


* The operator Cf{N) is regularized in the space and it changes the graph of a state by adding special 
loops at the nodes. Therefore, the dnal operator C^{N) acting is removing special loops at the nodes. 

® It was shown in [l6| that it is possible to obtain a self-adjoint curvature operator from the non-symmetric 
operator. 
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where 


0 


ei,ej 


{_^ijkJv,ejjJv,ej,k 

) {^ij'Jv,ei Jv,ej,k ^) ? 

2 tt 

Jv,ei ,j 

A/j 

_ y Jv,ei ,kJv,ej ^ky Jv,ej ,lJv,ej 


(IV.14) 

(IV.15) 


where k{v) 
free integer 


IS an averaging 
parameter [Ifil ]. 


coefficient that depends only on the valence of the node v, 
and V~^ is the “inverse volume” operator defined as 


A/j is a 


— + t^ll) ^V. 


(IV.16) 


The operator C^{N) is not graph changing and passes naturally to the space It maps 

its dense domain C "H^x to a subset of "H^x therefore preserves the gauge and partial 
diffeomorphism invariance. 


B. Quantum constraints algebra, symmetric constraint operator & physical states 
We can now introduce the non-symmetric scalar constraint operator 

C{N) := C^{N) + C\N) = ^ N{v){C^ + 4^) =: ^ N{v)C, . (IV.17) 

It is defined on a dense domain &{C{N)) C preserves "H^x- Since the classical 

scalar constraint functional is an observable, it is generally assumed that the quantum operator 
corresponding to it must be self-adjoint. However the operator in pV.17p is not symmetric and 
it was argued in [l^ that it is not necessary to have a self-adjoint constraint operator exactly 
because it is a constraint^®. Since we will be looking for the kernel of the scalar constraint 
operator, it may not be relevant to construct a self-adjoint operator as long as zero belongs 
to its spectrum. We will show below how we could introduce a symmetric constraint operator 
which is the first step toward defining a self-adjoint operator. 


1 . Quantum constraints algebra 

Let us for the moment assume that our constraint operator is C{N), then we can make a 
short calculation to check if this operator is anomaly free. The calculation goes as follows: given 
a state 'k.y G have 

[C{N),C{M)]^^= Ar(u)M(u')[4,4']'I^7 • (IV.18) 

v,v'£'y 

Because the regularization used to construct the operator is local with respect to each node, the 
commutator 


[ 4 , 4']^7 = 0 , y V ^v', 


(IV. 19) 


It was also shown in [l^ that a symmetric constraint operator may not be anomaly free. 
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hence 


[C{N),C{M)]^^ = Y,N{v)M{v)[C,,C,]% . (IV.20) 

vG'y 

In the space the commutator [Cy,Cy] also vanishes for the same reason in the case of 
Thiemann’s constraint operator, namely the two terms of the last commutator, when acting on 
a state in (before taking the limits of the regulators lIV.lip . produce two diffeomorphism 
equivalent states, therefore the commutator vanishes^^ on any state in "H^x 

[C{N),C{M)]=0 . (IV.21) 

When it comes to the algebra with respect to the other constraints, we already know that, 
on one hand, the operator C{N) preserves the SU{2) gauge invariance, on the other hand, 
a diffeomorphism constraint operator does not exist in this representation and the only thing 
we could check is whether it is covariant with respect to the action of diffeomorphisms. The 
calculation and the result is not different than in the case of Thiemann’s constraint operator 
and we find that indeed the operator C{N) is diffeomorphism covariant 

uy^[C{N)]Uf = c{f*N) , V / E Difr°°(s) . (iv.22) 

Therefore we conclude that the scalar constraint operator C{N) is anomaly free. 


2. Symmetric scalar constraint operator & physical states 

Concerning the question of defining a symmetric scalar constraint operator, it turns out that 
it is actually possible to introduce a symmetric operator using C{N) and its adjoint operator^^. 
The adjoint operator C'^(A^) is closed and also densely defined C ^(C^(N))), hence the 
operator C{N) is closable^^ and (C'^(A^))^ = C{N). Therefore in the rest of the article we 
consider the closure of C{N) and C^{N) as being the non-symmetric scalar constraint operators 
at our disposal. 

The operator C^N) could be by itself considered as a quantization of the classical scalar 
constraint functional IIV.II and it could stand as the quantum scalar constraint operator in the 
theory on the same footing as the operator C{N). If the implementation of the scalar constraint 
is appropriate, then in the semi-classical limit of the theory the expectation values of the operator 
and its adjoint should coincide, up to small quantum corrections. Hence, both operators are 
equally good candidates for the scalar constraint operator in the theory. Notice that C^{N) is 
also anomaly free, i.e. it preserves SU{2) gauge invariance and we have 

[C^{N),C^M)] = 0 , Uj^[C^iN)]Uf = C^{f*N) , V / E Diff“(S) . (IV.23) 

In order to construct a symmetric scalar constraint operator Cgyrai^), we suggest to define 


The commutator vanishes with respect to URST (topology) [rI. [lll|. 

Definition: Let T be a densely defined linear operator on a Hilbert space Jif. Let ^(T^) be the set of G 
for which there is an G with 

(T’b, fp) = (b: h) for all "b E t^{T) 

For each such p G we define = rj. The operator is called the adjoint of T. 

We keep the same notation for C{N) and its closure. 
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it as a combination of C{N) and C^{N). The simplest example is 

aym{N) := ^{C{N) + C^N)) , ^(Csy^iN)) = &{C{N)) n 9{C\N)) . (IV.24) 

It is obvious that this operator is closable, densely dehned and anomaly free. The question 
of existence of self-adjoint extensions is still open. However it is a strongly eligible candidate 
for the scalar constraint operator in the theory. The structure of its kernel, equivalently the 
solutions to this constraint in the space can to some extent be described quite easily. The 

properties we know so far of the kernel elements of C{N) and C'^(A^) can be summarized as 
follows: 

• every state that is in the kernel of the volume operator V and has coplanar edges at all 
the veritices of its graph, is in the kernels of C{N) and C'^(A^); 

• the set of states of non-zero volume^^ in the kernel of C{N) contains an infinite number 
of states that have the form of finite linear combinations of spin network states^®; 

• states of non-zero volume that are in the kernel of C'(A^)^ have the form of infinite linear 
combinations of spin network states; 

• states of non-zero volume with graphs that do not contain special loops are neither in the 
kernel of C{N) nor the kernel of C'^(A^). 

those properties, we can deduce that the kernel of C'sym(fV) has the following structure: 

every state that is in the kernel of the volume operator V and has coplanar edges at all 
the veritices of its graph, is in the kernel of C'sym(iV); 

states of non-zero volume that are in the kernel of C'sym(fV) have the form of infinite linear 
combinations of spin network states: 

• states of non-zero volume with graphs that do not contain special loops are not in the 
kernel of Cgymi^)- 

Having a scalar constraint operator, as we mentioned before, the construction of physical 
states is achieved via averaging of the elements of its kernel, subset of with respect to the 
rest of diffeomorphisms in Diff(Il)/Diff(S)vert(7)- 


With 


V. COMMENTS & OUTLOOKS 


In this article, we presented a concrete implementation of the scalar constraint operator in 
loop quantum gravity. The construction of the Euclidean part of the constraint operator uses 
a regularization based on the assignment of ’’special” loops 0 , while for the Lorentzian part 


of the constraint we use the curvature operator of 1^. The resulting non-symmetric operator 
C{N) is densely defined on the Hilbert space of partially diffeomorphism invariant states 
introduced in 1^. The operator C{N) is SU{2) gauge invariant and diffeomorphism covariant, 
it preserves the space and its algebra is anomaly-free. 

Thanks to the properties of the special loops, the adjoint C'^(iV) is a densely defined operator 
on has the same properties as C{N). It also allows to construct symmetric constraint 

operators, Csym{N), as combinations of the operators C{N) and C^N). The operators C{N), 


By a state of non-zero volume we mean any state which is not in the kernel of the volume operator. 
Simple examples of such states can be straightforwardly derived. 
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C'^(A^) and Csym{N) are all equally suitable candidates for the scalar constraint operator in loop 
quantum gravity. In each case, the general structure of the kernel of the constraint operator is 
known on a qualitative level, as outlined in section IIV B 21 

The regularization proposed in this article could also be applied in order to define a Master 
constraint operator, corresponding to the classical Master constraint functional introduced by 
Thiemann [1^ as a way of reformulating the singular scalar constraints C(x) of equation pi.3p . 
Carrying out the construction, one would obtain a densely defined operator on which is 

symmetric, gauge and diffeomorphism invariant, and anomaly-free. However, further work is 
needed in order to investigate the structure of the kernel of this operator. 

The freedom of choice between different eligible scalar constraint operator should be regarded 
as a quantization ambiguity that can be fixed only through a semi-classical analysis of the 
dynamics in the theory. Therefore the next step of our program is the challenging task of 
constructing, or at least approximating, semi-classical states in the theory. The example of the 
operator C{N) is encouraging in this direction since its kernel is more tractable with respect to 
the spin network basis. 
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